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Abstract 

It is demonstrated that almost any S-matrix of quantum field theory in curved 
spaces posses an infinite set of complex poles (or branch cuts). These poles can 
be transformed into complex eigenvalues, the corresponding eigenvectors being 
Gamow vectors. All this formalism, which is heuristic in ordinary Hilbert space, 
becomes a rigorous one within the framework of a properly chosen rigged Hilbert 
space. Then complex eigenvalues produce damping or growing factors. It is known 
that the growth of entropy, decoherence, and the appearance of correlations, occur 
in the universe evolution, but only under a restricted set of initial conditions. It is 
proved that the damping factors allow to enlarge this set up to almost any initial 
conditions. 
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1 Introduction 



For many years unstable quantum states were represented by Gamow vectors Q, i.e., 
eigenvectors corresponding to complex eigenvalues of the hamiltonian. But since the 
hamiltonian is self adjoint, if we use a Hilbert space as state space, eigenvalues must be 
real. For this reason, Gamow vectors were deshonorably excluded from ordinary quan- 
tum mechanics and they were considered just as useful (but not rigorous) analogies or 
approximations. Nevertheless some years ago it was proved that Gamow vectors belong 
to an extension of Hilbert space, namely a rigged Hilbert space with a nuclear subspace 
based in Hardy class functions (see e.g. 0, [| |J and bibliography therein). Since then, 
Gamow vectores are legal citizens of an extended version of quantum mechanics, where 
complex eigenvalues greately help the computation of survival probabilities, life times, 
Liapunov variables, the evolution toward equilibrium, etc. [H, ||. These eigenvalues also 
allow the introduction of more refined physical concepts, e. g. the thermodynamical 
arrow of time can be defined in a way which is free of the usual criticisms (Lochshmidt 
objection, coarse-graining ambiguities, non-systematic approximations, etc.). Precisely, 
studying this arrow of time we can see that all the problem of time asymmetry essentialy 
has a cosmological origin |5|, |Bj . Therefore it is natural to demand if unstable quantum 
states, considered as vectors of a rigged Hilbert space, can be used in quantum cos- 
mology. This is, in fact, the case and the first examples are Refs J7[ and j8| where a 
simplified (toy model) version of the universe is studied as a Friedrichs model [4]. 

In this paper we will present a complete (not only a toy model) semiclassical model of 
the universe following the line of Refs. || and jnj and we will show how the presence of 
unstable quantum states enlarges the set of cases where we can prove that the decoher- 
ence phenomena appears. Also correlations appear (for unstable states) explaining the 
outcome of a classical universe. Essentially, in this paper, we study, just one example, 
but we will also especulate about eventual generalizations. 

But first, let us briefly recall the guiding lines of the extension from Hilbert space 
to rigged Hilbert space in usual quantum mechanics. The traditional set of states of a 
quantum system is a Hilbert space, which leads, as it is well known, to time-reversibility. 
It is precisely this property which changes drastically with the extension of the Hilbert 
space TC to a rigged Hilbert space. This extension corresponds essentialy to the transition 
from a space of square integrable functions to a space of distributions. This procedure is 
not unique and different distribution spaces can be defined which are based on different 
test function spaces. If we choose as the test function space generated by the 
eigenfunctions of the energy uj which are analytic in the lower complex halfplane, when 
the real variable uo is promoted to a complex variable z (precisely Hardy class functions), 
we obtain the dual space $*, which is the required extension of the Hilbert space Ti,. 
The corresponding Gel 'f and triplet is then 

$_cftc$*. (l.i) 

If the same procedure is performed in the upper complex plane, the resulting triplet 
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reads 



$+CftC$*. (1.2) 

As we will see the first of these choices, hence the space <&*, corresponds to unstable 
dacaying states while the second one, namely $* , corresponds to unstable growing states. 
In fact, the complex poles of the S-matrix are related, as it is well known, with unstable 
physical states. These poles can then be transformed into complex eigenvalues z n of the 
Hamiltonian using standard methods [3]. The essence of the rigged Hilbert space rather 
than the Hilbert space, as the framework for the quantum states of the system, is clearly 
exhibited precisely at this stage: the eigenvalues z n of a Hermitian operator are not real 
anymore in this extended space |], |TT|. If Im z n > then a growing prefactor appears in 
the time evolution of the corresponding eigenvector \n + ), giving rise to a growing state 
belonging to the rigged Hilbert space <i>^ . On the contrary, if Im z n < the prefactor is 
a decaying one, the corresponding state is decaying and belongs to another rigged 
Hilbert space, $*. Finally, Im z n = corresponds to an ordinary stable state belonging 
to the ordinary Hilbert space 7i = D (more general models contain both, growing 
and decaying states [5,6]). 

If K is the Wigner time-reversal operator we have 

K : $* -»• ■ K : -»• (1.3) 

since decaying states are transformed into growing states (or vice-versa) by time-invers- 
ion. Then the choice of (or $*) as our space of quantum states implies that K is 
not defined inside (or $+), so that irreversibility naturally appears and therefore the 
arrow of time also appears in the quantum regime. 

It follows that the choice between $_ or $ + is irrelevant, since these two objetcs 
are identical (namely one can be obtained from the other by a mathematical symmetry 
transformation), and therefore the universes, that we will obtain with one choice or the 
other, are also identical and not distinguishable. Only the names past and future or 
decaying and growing will change but physics is the same and, e.g., we will always have 
equilibrium toward the future. 

Let us summarize the organization of this paper. Section 2 introduce the model that 
we will study. In Section 3 it will be demostrated that the S-matrix of the model has an 
infinite set of complex poles (or branch cuts), and how these poles are transformed in 
complex eigenvalues that originate, in turn, damping factors of unstable decaying states. 
In Section 4 it is shown how these damping factors enlarge the set of initial conditions 
where decoherence occurs. In Section 5 we will prove that there is correlation in all 
unstable states. Finally, we briefly state our conclusions in Section 6. Two Appendixes 
complement this work. 
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2 The Model 



Let us consider the model of Sec. 3 of Ref. [10] where a Robert son- Walker metric is 
studied (that we will mainly consider in the flat case), with a total action S = S g + Sf, 
being S g the gravitational action and Sf the matter action (the usual action of a spinless 
massive field $). The gravitational action is given by 



S g = M 2 J 'd V [~a 2 -V(a)\, (2.1) 



where M is the Planck mass, r\ is the conformal time, a is the Robertson- Walker 
scale factor, a = da/drj, and V(a) is the potential function that arises from the spatial 
curvature, a possible cosmological constant and eventually a classical matter field. As 
this last field is arbitrary, for the sake of simplicity, let us study the case where the 
classical matter field is such that V(a) = B 2 / 2(1 — A 2 /a 2 ) where A and B are arbitrary 
constants. 

This case is the simplest of all, but we believe that the main features that we will find 
will also be presented in more general cases, as we will argue bellow. The role played 
by the classical field is completely natural, in the context of this paper. In fact, we 
will essentially work using some results of quantum field theory in curved space-time, 
where the geometry of space-time is fixed "a priori" (namely there is no back-reaction). 
The classical field is, precisely, the agency that do this job, fixing a class of possible 
geometries (but the properties that we will find will be the same for almost all classical 
field). 

The Wheeler-DeWitt equation for our Robertson- Walker model is 

'^d a * + M 2 V(a) - ± jf [d^ - C^f)]*(a, $) = 0. (2.2) 
Thus after making the WKB ansatz, the Hamilton- Jacobi equation appears as [10[ 
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where S is the principal Jacobi function. Thus the (semi) classical time parameter or 
WKB time rj is given by 



(2.4) 



d dS d 
drj da da 

Then in our simplified model we have the following class of geometries, in terms of 
this conformal time drj = a _1 dt [|i~2 



a = ±(A 2 + B 2 r] 2 )2 + C (2.5) 

where C is an arbitrary constant. Using different values for this constant and different 
choices for the ± sign we obtain different classical solutions (in a more general case many 
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constants would be necessary). Going now to Ref. [|13j (eq. (3.113)) we can see that 
the semiclassical (or quantum field theory in curved space-time) problem is solved for 
all four dimensional, spatially flat, cosmological models with scale factor 



C(rj) = a 2 = A 2 + B 2 rf - oo < r] < oo (2.6) 
where A and B are constants. Then if we consider a massive, conformally coupled scalar 



field, the energy function VL 2 k reads 



n 2 k = m 2 a 2 + k 2 = k 2 + m 2 (A 2 + B 2 V 2 ) (2.7) 

where m is the mass of the quantum matter field and k 2 = \k\ 2 , where k/a is the 
linear momentum of this field, in the case of flat space Robertson- Walker universe (or a 
function of this momentum in the two other cases, namely open and close, being k is a 
discrete variable in the close case). Then (2.7) coincides with the last equation of page 
70 of Ref. |J. 



If we ideally consider the evolution of the universe from rj — > — oo to rj — ► +oo (even 
if really we would like to have only an expanding universe and therefore rj > 0, we 
will discuss this issue below) and we define the corresponding adiabatic vacua |0,m) 



for rj — > — oo and |0, out) for rj — > +oo the Bogolyubov coefficients are (Ref. |L3|1 , eq. 
(3.124)) 

a kj = r[ i (1 _ iAfc)] S* = ( 2 -8) 

7T 

(3 kj = -iexp(--\ k )5 kj = f3 k 8 kj , (2.9) 

where X k = k 2 /Bm + A 2 m/B and 5 k j is the Kronecker S, for the discrete case, and the 
Dirac 5, for the continuous one. 

Let us comment now on the choice of the vacua since really we would like to study 
only the evolution 77 > 0. The |0,out) vacuum is the adiabatic physical vacuum for 
a — > +00, where the classical regime must naturally appear, therefore it is a completely 
reasonable vacuum. Let us suppose that the vacuum at r] = is just \0,in). This is, of 
course, a completely arbitrary choice, that will be discussed in the next section, where 
we will introduce a general vacuum at rj — that we will call |0, 0). Anyhow with this 
arbitrary choice (2.8) and (2.9) are corrects. 



The poles of the S-matrix and the unstable quan- 
tum states 



From (3.46) and (3.47) of Ref. [13|, or more generally from Sec. 2 of Ref. it can be 



seen that there is a pole in the S-matrix (between the "in" and the "out" Fock spaces) 
where the function Aj-j = —iY, k PkjCtik has a pole, namely where a k j = ( or (3 k j has a 
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pole). Using (2.8) it must be = or, which is the same thing, that T[l/2(1 — iXk)] 
would have a pole. T(z) has a pole if z = —n (n= 0, 1, 2,..., see, e.g., [jnj or [|H|) (no 
poles are produced by the /3's given by (2.9)). Therefore S has a pole if 




2i(n+~)\, 



(3.1) 



and the squared energy, for each pole, reads 



n 2 k = m 2 a 2 + mB[ 



mA 2 



2i(n+-)}. 



(3.2) 



B 



We will call this energy Vt^ simply Q n . Thus, we have an infinite set of unstable 
states with mean life 



Let us observe that the energy and mean life are a-dependent. Therefore we have 
two possibilities: 

1) either we can consider that the in and out states corresponds to a » 1, where 
these mean lifes are big but still finite, or 

2) we transform all the equations to the non-rescaled case, where the physical real 
values are the physical time t = J adrj, the physical energy Qk/a and the physical 
momentun k/a. 

We follow the first alternative, and sketch the second one in the Appendix A. 

Therefore the universe evolution creates unstable particles as well as stable ones. 
Using the standard method explained in Refs. [Q] and |J we can promote these unstable 
states to vectors of an adequate rigged Hilbert space and build a basis of this space with 
stable modes with real energies Qk plus unstable modes with complex "energy" Q n given 
by (3.2) (in the open case this procedure is direct, since we have a continuous spectrum 
to begin with, in the close case we must use assumption 3 of Ref. J7J). 

This would be the state of affairs if we use the (quite arbitry) vacuum \0,in) of 
section 2. In this case we have found an infinite discrete set of unstable states. What 
happens if we use a generic (i.e., almost any) vacuum 1 , 0) at rj = 0? This generic 
vacuum will be related to |0, in) by some Bogolyubov coefficients dkj, Pkj- Then the a 
coefficient relating |0,0) to \0,out) reads 



The poles are now located where this alpha vanishes. The roots in k of the corre- 
sponding equation, = can be found only if we fix the arbitrary coefficients a^, fiik- 
Of course if these coefficients are fixed in a very particular way the equation will have 
no roots. But if the functions a:^, (5^ are fixed in a generic (i.e. in almost any) way, the 



2i{m 2 (a 2 - A 2 ) + [m 4 (a 2 - A 2 ) + AB 2 m 2 



jn + \) 2 ]^ 



(3.3) 
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equation will have a set of complex roots, that correspond to unstable particles created 
by the universe evolution. This statement is equivalent to claim that a generic S-matrix, 
for our problem, has infinite numbers of poles or cuts. Let us say an infinite set of poles 
to precise the ideas (cuts will be studied in Appendix A). Even if we have not, by now, a 
rigorous mathematical proof of this theorem, we think that we can sketch a reasonable 
convincing demonstration 

Let {\n, in)}, (|0,0)}, and {\m, out)} be the basis of the Fock spaces corresponding 
to vacua |0, in), |0,0), and |0, out). The in-out S-matrix, with an infinite set of poles, 
reads 

Sum = (n,in\m,out) = *^2(n,in\l,0)(l,0\m,out) (3.5) 

i 

where the states |/,0) are a complete set. From some (infinite) values of n and m we 
know that S nm has poles. Let us consider one of these values, then the l.h.s. of (3.4) 
has also a pole, and theferore one of its terms has a pole. Then, either: 

i) one of the factors inside the summatory of (3.4) has a finite number of poles and 
the other one an infinite number of poles, or 

ii) both factors have an infinite number of poles. 

But (i) must be excluded since, in this case, time evolution (— oo < r\ < 0) would 
be qualitatively different to evolution (0 < r\ < +oo) and this fact would break the 
time symmetry, which is impossible since evolution equations, time evolution of a, and 
boundary conditions are time symmetric with respect i] = 0. Then, the — out matrix 
(l,0\m, out), corresponding to evolution i] > 0, has an infinite number of poles. 

We give an alternative demonstration in Appendix B, which is valid for every evolu- 
tion and every spatial geometry. 

Even if a rigorous proof of these facts would be welcomed we believe that the reason- 
ings above, and the ones in Appendix B, are quite convincing. Essentially the periodical 
nature of ot^ is inherited from its definition; = (ui,Uk) (eq. (3.36), Ref. [13]) where 
Ui and Uk are two different negative frequency solutions of the corresponding Klein- 
Gordon equation. As in flat space-time, these solutions are functions like exp(— ikt), 
they somehow must keep the periodicity in k in curved space-time. The a coeficients 
always have a periodic behaviour in the complex plane as it is shown in equations (3.91), 
(3.124), (4.60), (4.61), (4.95), (5.41), (5.110), and (5.111) of Ref. [13]. Therefore the 
S-matrix has an infinite and discrete set of complex eigenvalues for almost any initial 
condition |0, 0). 

Now that we know that the 0-out S-matrix has an infinite set of complex poles, we 
can find the complex eigenvalues [2,3]. 

As the { \k, out)} basis is complete we have 

/ | k, out){k, out \dk = 1, (3.6) 

J k 

where \k,out)<E H and the integral means that we must integrate over the continuous 
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spectrum of energies and other quantum numbers. Using the standard techniques of 
Ref. [2,3] we can transform the last equation in 



^2 l n J out— ) (n, out + | + / \k, out— )(k, out + \dk — 1, (3.7) 

where \n, out— ), \k, out— )e and the first summatory corresponds to the discrete 
unstable modes and the integral to the stable continuous ones. 

We choose a $_ test function based in Hardy functions for below, all the poles will 
have negative imaginary part and all the unstable states are decaying ones, and they 
will belong to (as we already know we can also make a symmetric choice). 

According to Wheeler-DeWitt equation (2.2), the field hamiltonian reads 

h = \l k (~ d k + n l$) dk = I Mk4a k dk, (3.8) 

where at and aj. are the usual creation and anihilation operators. From now, we will 
always refer to the out case with a » 1 and h, Qk, &k, and a\ will be h out ,fll ut , a>l ut , and 
a}^ ut . There are new creation and anihilation operators for the discrete spectrum: a° ut , 
o*out anc j f or con tinous ones a^ ut , a>l° ut (the definition of * is given in Appendix B). 
Vectors \n, out—) will be created by the repeated action of a* ou ' on |0, out), and vectors 
\k, out— ) will be created by al° ut analogously. 
Therefore h out now reads 

h out = J2K° ut < t + I Sl k al out a° k ut dk, (3.9) 

and we will have 



h out \n, out) = Q n n\n, out), (3.10) 

so the evolution of \n, out) e $^ has a damping prefactor exp(— ni]/T n ) since fl n has an 
imaginary commponent. 

Thus, as we now have damping factors exp(—nr)/r n ) in the evolution equations, it will 
be very easy to find Lyapunov variables, and in particular a growing entropy for almost 
any initial condition as in Refs. || [|. This result must be compared with the one of Ref. 
T7| where a Luapunov variable was found for the universe evolution using the standard 



methods |18|, [L9[ based in an arbitrary coarse-graining and a particular (generalized 
molecular chaos) initial condition. The new result is much more satisfactory than the 
old one, since now we have a growing entropy for almost any initial conditions solving 
Lochsmidt criticisms. 



4 Decoherence 



Decoherence is a dissipative process, and we know [20| that it is closely related to 



another dissipative phenomenon, that is particle creation from the gravitational field 
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Particle creation has been studied in the quantum field theory in curved 



to the semiclassical approximation to study the decoherence phenomenon. 

Decoherence naturally appears in systems where the hamiltonian has complex eigen- 
values, as it is proved in Ref. 



25 



and [|K] 



Let us consider the formalism developed in Refs. || 
to see that, this is also the case, in the system we are studying. We labelled 
the three-geometry with the scalar factor a (the indices a and (3 simbolizes the choice of 
the sign and constant in (2.5)), and $at is the mode N of the matter field; precisely, we 
have used n for the discrete unstable states and k for the continous stable states; when 
we will be refering to both kinds of states, we will call the index N. 

The WKB solution of the Wheeler-DeWitt equation reads (|| eq. (2.8)) 

*(a, = exp[^M^(a)] X (a, [$ N ]), (4.1) 

where S is the principal Jacobi function of (2.3) and x( a > [*&n]) can always be written 
as 



X (a, [Q N ]) = Y[ X n(v,®n). 

N 

We can obtain Xn(Vj ®n) by a Gaussian approximation [|1(J 
Xn(v, ®n) = A N (rj) exp[ia(r]) - B N (rj)$ 2 N ]. 



(4.2) 



(4.3) 



Functions A N {ri) and ocn(v) are rea l while B N {rf) is complex, precisely B N (rj) 
Bnr(v) +iB NI (rj) and can be obtained solving the system flO|j : 



7T 4 



(2B NR ( V ))- 



(4.4) 



ctN(v) = -B N r{v), (4-5) 



B N ( V ) = -2iB 2 N (r ] ) + i -n%( V ). (4.6) 



2 

From [10] or [20] we can learn the conditions for the ocurrence of decoherence if we 



use only the real fljy, namely the ones of the continuous spectrum. In Ref. |26[ the 
computation is made only through a linear approximation of B(a) as a function of a. In 
Ref. [20] decoherence takes place only in the case where the Bogolyubov coefficients f3 n 
are small and imaginary. In Ref. [10] decoherence takes place unless the environment is 
very ordered and fine tuned. Thus we cannot say that there is decoherence for almost 
any initial condition if the fl^ are all real, like in these works. Let us see, what happens 
in our model if we use a basis with infinite complex modes Q n as well as real modes fife. 

From the wave function (4.1), and after the integration on modes of the scalar field 
(considered here as the "environment"), we obtain the following reduced density matrix 



8 



p r (a, a') = exp[—iMS a (a) + iMS a (a')]p aa (a, a') 
+ exp[-iMS a (a) + iMSp(a')}p a/3 (a, a') 
+ exp[-iMSp(a) +iMS a (a')}p f3a (a,a') 

+ exp[-iMS,3(a) + iMS p (a')]p^(a, a'), (4.7) 

where, as we have said, a and (5 symbolize two different classical solutions, namely two 
different choices of the sign ± and the constant C of (2.5), and 

pf(a,a') = J{pf N (a,a') =\{ f d<S> N X a N *(v,®N)xM,®N). (4.8) 

N N J 

From (3.20) of Ref. [10], it is 

B N = -i^, (4.9) 

where g n is the wave function that represents the quantum state of the universe being 
also the solution of the differential equation 

-g N + n 2 N g N = 0, (4.10) 

Qjy can be the complex energy Q n in our treatment. From (2.6), (2.7), and (3.1) we 
know that if the initial sate is |0, in), the complex energies are 

n*=m 2 a 2 -rn[2iB(n + -)+mA 2 ]. (4.11) 

In the more general case we use an arbitrary initial state |0,0), instead of |0, in). 
From the discussion of Sec. 3, we know that, in a generic case, an infinite set of complex 

poles does exist. Then we must change (3.1) by k 2 = (n = 0, 1, 2, ) where this are 

the points where the infinite poles are located in complex plane k 2 ; (3.2) now reads 

n 2 n = m 2 a 2 + k 2 n . (4.12) 

Let us now see that decoherence takes place if there is an infinite set of complex 
modes (even in a more general case than the one of the time evolution fixed by (2.3), 
the only one we have studied in great detail above, if we use the theorem of Appendix 
B). 

Let us consider the asympthotic (or adiabatic) expansion of function g N when a — > 
+oo in the basis of the out modes. g N is the wave function that represent the state of 
the universe, corresponding to the arbitrary initial state |0,0), and its expansion reads 
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9n = /S- ex PH / ^Ndr]} + -Q^=exp[i [ Vt N drf[ : (4.13) 

where Pjy and Qat are arbitrary coefficients showing that |0, 0) is really arbitrary. 

It is obvious that if all the are real, like in the case of the f^> (4-13) will have 
an oscilatory nature, as well as its derivative. This will also be the behaviour of in 
(4.9). Therefore the limit of B^ when 77 — > +00 will be not well defined even if B^ itself 
is bounded. 

But if Qn is complex the second term of (4.13) will have a damping factor and the 
first a growing one. In fact, the complex extension of eq. (4.13) (with N = k) reads 

9n = ^- exp[-t| n n dr]] + -^=exp[i^ ^ n drj], (4.14) 
Therefore when 77 —>■ +00 we have 

5n « = V- (4.15) 

Then we have two cases: 

i) = fifc G 1Z + for the real factors. Then we see that when rj — > +00, the r.h.s. of 
(4.8) is an oscillatory function with no limit in general. We only have a good limit for 
some particular initial conditions (as Qn = or P/v = [§, |H], [2(| |26fl). 



ii) f2jv = = E n — f^ 1 € C for the complex factors. If we choose the lower Hardy 
class space $_ to define our rigged Hilbert space we will have a positive imaginary part, 
and there will be a growing factor in the first term of (4.13) and a damping factor in 
the second one. In this case, for a — ► +00, we have a definite limit: B n = 1/2 Q n . 

So we can say nothing about the limit of the real factors (and therefore nothing in 
general for the product of these real factors) while the complex factors have definite 
limits for every initial conditions, namely for every |0,0) state. 

Therefore let us compute the p^f for the complex factor, since these are the only 
quantities whose limits we know for sure. So let us compute these matrix elements using 
eq. (2.29) of Ref . [9] or (2.24) of [10], namely 

P ™ {a ' a)= \ [B- n {r,,a) + B n (TfMV ^ «<n> a ) + ta »to>M> (4 ' 16) 

where a and (3 mean that the B refers to these classical solutions. p^(a,a') can be 
obtained using (4.8). Let us first compute 

log\pf (a,a')\ = £ log\p*(a, a% (4.17) 

n 

Now it can be proved that if TmB n »s 2m^Q n 7^ 0, 
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\p£{a,a!)\<l. (4.18) 
In fact, calling B^(r], a) = z = x + iy and B n (r)', (3') = ( = £ + iy, we can compute 



4< | < 4x£ 



{x + C) 2 + {y + r])] 2 - \ X + S\2 ^' ^ 19 ) 

since from \x + £| 2 > it follows that AxS, < \x + £| 2 . Then all terms of the r.h.s. of 
(4.17) are negative if Im B n ^ 0. 

Now from (4.11) or (4.12) we can see that the B n (r], a) ~ l/2Q n (i], a), corresponding 
to the discrete complex modes, have all almost the same asympthotic value when a — > 
+oo, therefore all the terms of the r.h.s. of (4.17) have almost the same asymthotic value. 
As they are all negative and almost equal, the summatory of (4.17) has an asynthotic 
value — oo and therefore p° 'P(a, a') of (4.8) vanishes only considering the discrete complex 
modes factors. 

Then we have decoherence if B*(rj,a) ^ B n (r}',/3) namely if f2*(?y, a) ^ £l n (r]',f3), or 
using (4.12), if (for an infinite set of n) we have 

m 2 [±(A 2 + B 2 rf)^ + C a ] 2 + m 2 [±(A 2 + B 2 V ' 2 )^ + C p ] 2 . (4.20) 
So we necessarily have decoherence: 

i) for different classical solutions, i.e. for different constants C a ^ Cp, or different ± 
signs, even if the time is the same 77 = i]'. 

ii) for the same classical solutions (C a = Cp, and same ± sign) if the times 77 and rj' 
are different. 

Now we can discuss the choice of either <3>_ or $ + for the test function spaces. $_ 
produces the space $^ of the decaying states, while $ + will produce the space $^ or 
growing states. Of course one choice becomes the other if we change the ± sign in (2.5) 
and also in all other equations that are a consequence of (2.5). To choose $_ or $ + 
corresponds to the choice of the arrow of time in the quantum regime as explained in 
the introducction. Choosing the Gel'fand triplet $_ C T~i C is equivalent to say that 
the unstable created particles produced by the universe expansion will decay. This is 
the motivation of the choice of the arrow of time in the quantum regime. 

Nevertheless we can conceive more complex models where growing and decaying 
particles could coexist at the same time. Then, we can adopt a more conservative 
attitude. We can consider that the space $ + ©$_ as the space of the test functions with 
the correponding space <3>^ © $^ where there are mixed decaying and growing states. 
This choice is possible in the quantum regime, and therefore there will be no arrow of 
time in this regime. But in the classical limit, solution in space $^ will decohere with the 
solutions in $^ since they correspond to different choices of the ± sign [26]. Therefore 
for a classical universe either we have a state in $^ or $^ (which, on the other hand, 
it is an irrelevant choice). Under this perspective, we will not have an arrow of time in 
the quantum regime but this arrow will appear naturally in the classical regime. 
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5 Correlation 



From Ref. [10] we can also learn the conditions for the existence of correlation. But, for 
the same reason used in Sec. 4, we cannot say that there is correlation for almost any 
initial condition. This correlation depends on the initial conditions, and it can be easily 
obtained from system (4.4)-(4.6). In fact for certain initial condition, if the Q% are real, 
it turns out that B^r — and all the conditions for the existence of correlations of Ref. 
[10] are not fulfilled and therefore there is no correlation. Precisely, if B^r = and 
B N = iB NI and all energies fl^ are real (it would be better to call it flk), (4.6) reads 



B 



NI 



2B 2 NI 



2 



N • 



(5.1) 



where all the variables are real. Therefore if B^r = at a = 0, B^r = at every time 
and there is no correlation. Thus correlation, as decoherence, depends crucially on the 
initial condition and there is no correlation for the above initial condition. Correlation 
takes place inside each classical solution and it therefore can be computed using the 
Wigner function, associated with PrniViV') EH EI 



[n)(a,P) = j dAexp(-2iPA)/C(a-— ,a + 



(5.2) 



a ± -tt, and P is the canonical momentum. 



where a, a' 

Nothing new can be said about the real continuous modes, all was already said in 
Ref. [10]. We must only study the complex discrete unstable modes. This is nevertheless 
important since, most likely, the universe is in an unstable mode or more generally in a 
linear combination of unstables modes (see Ref. [7] and |28| , |29fl , where the universe is 
in a "tunneling" unstable state, i.e. a typical Gamow vector). 

Then we can repeat the reasonings of [10] from (2.24) to (2.28) and, with the same 
assumptions we will arrive to this last equation, that now reads 



7T 



F« a (n)(T,P)^C 2 (T)J-exp 



(P-MS + a- 



4B„ 



(7 Z 



In the case of the our unstable states we have, for a{rj) — > +oo 



(5.3) 



a 



B n R — T=- 



V2 



m z B 2 7] 2 + (m 4 5V + ^n A B\n + -) 



(5.4) 



B nI = 2V2- 



m B 77 



1 

2l2 



m 2 B 2 Tj 2 + ^ m 4_BV + 4m 2 B 2 (n + |) 2 J J + Am?B 2 (n + \) 2 
and the inverse of the correlation width of the reduced density matrix is 



(5.5) 
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a 2 = 



m^B 4 ^ 2 \m 2 B 2 rf + (m 4 B 4 r] 4 + Am 2 B 2 (n + \)) 



m 2 B 2 r] 2 + {jrfB^rf + Am 2 B 2 (n + + Am 2 B 2 (n + ±) 



(5.6) 



if 77 is big: \B R \ > \B 



We can see that when 77 



-00, a 



0, and there is a good correlation and the 



Wigner function is a Gaussian function, of width a, peaked about 

B n i 



P = MS-a + 



4B. 



nR 



(5.7) 



where the first term of the r.h.s. gives the classical result and the last two are the 
quantum correlation to the classical trajectory. 

On the other hand, in this state we can predict strong correlations between coordi- 
nates and momenta [10], because 



MS-d + 



B 



nl 



AB. 



nR 



» a 



(5.8) 



(in the preceeding equation we can see that, in the limit of large 77, the l.h.s. is propor- 
tional to i] 2 , while a 2 ~ t^ 1 ). 

For the generic initial condition |0, 0) we can use (4.12) and we will reach to the same 
conclusion. Therefore there is a perfect correlation for the unstable states of our model. 
Decoherence and correlation will produce the outcome of an asymthotic classical regime 
in the far future. 



6 Conclusions 

We have demonstrated that, the S- matrix of almost any quantum field theory in curved 
spaces model, has an infinite set of poles (or cuts). The presence of this singularities 
produces the appearance of unstable states (with complex eigenvalues) in the universe 
evolution. The corresponding eigenvectors are Gamow vectors and produce exponen- 
tially decaying terms as in the Friedrichs model of resonances. But the best feature of 
these decaying terms is that they simplify and clarify calculations. 

E. Calzetta and F. Mazzitelli [20] have demonstrated that, under suitable conditions, 
the expansion of the universe leads to decoherence if this expansion produces particle 
creation as well. Our unstable states enlarge the set of initial conditions where decoher- 
ence occurs. In fact, the damping factors (related to the imaginary part of S-matrix's 
poles), allow that the interference elements of the reduced density matrix, dissapear for 
almost any initial conditions. 

Following the reasonings of Ref. [10], we also demostrate that the unstable states 
satisfy the correlation conditions, which, with the decoherence phenomenon, are the 
origin of the semiclassical Einstein equations. 
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For simplicity, we assume (as usual) that the state of the environment can be de- 
scribed by a Gaussian wave function (eq. (4.3)). This is indeed a restricted class of 
states [10], but general states could also be implemented in our formalism. The arbi- 
trary election of the coefficients Pjv and Qn, in eq. (4.13), shows that the set of initial 
conditions is really arbitrary. 

Finally, we can say that the existence of unstable states in the universe evolution 
(coming from singularities in the Riemann second sheet of the analytical extension of 
the S-matrix) can help us to understand the quantum to classical transition and other 
dissipative aspects of the universe evolution. 
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Appendix A 



Equation (2.7) reads: 

Ql = k 2 + m 2 a 2 , 
so the relation between the physical energy and momentum is 



^i-Y = m > + ( k -)\ 

From (3.1) we know that we will have a resonance if 

mB 



= to + 

a J 



mA 2 , 1 . 
- r -2,(n+-) 



(A.l) 



(A.2) 



(A.3) 



Since a — > +oo the only way to have resonance for energies different from to, is that 
at the same time n — > +oo, so let us define 



N(n,a) 



n 

^2 ' 



(A.4) 



when both a — > +oo and n — > +oo. If n grows in the proper way, N(n, a) remains finite 
and we have 
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M = m 2 - 2imBN. (A.5) 



Then 



a 



1 

V2 



m 2 + Vm 4 + Am 2 B 2 N 2 



2mBN 



V2 



m 



2 + Vm 4 + Am 2 B 2 N 2 



(A.6) 



So the physical decaying time is 

V2a 2 



2mBn 



m 2 + \ m 4 + Am 2 B 



(A.7) 



where a >> 1 and therefore also n» 1 to get a finite result. 

Let us now compute the difference between the square of two subsequent physical 
energies 



n 



n+l 



Q ^ 2 



2mB 



(A.i 



When a — > +oo this difference vanishes and therefore, most likely, we have a cut 
singularity located at the physical energies 



a 



(m 2 - 2imBN)z. 



(A.9) 



It is easier to deal with a set of infinite poles that to work with cuts. Therefore we 
postpone, further discussion of cuts for future papers. 



Appendix B 

Let us demostrate that there are an infinite set of poles or cuts in the general case. 

Ordinary anihilation and creation operators are related by the usual canonical com- 
mutation relations 

[a u ,al,] = 5(uj - lo'). (B.l) 

If we promote the real variable uj to a complex variable z, it is demonstrated in Ref. 
0] that these relations must be subtituted by 

[a g ,a*,]=5(z-J), (B.2) 
where z, z' G T, T is a contour in the complex plane (for details see Ref. |J) and 

a* = Ka\K\ (B.3) 
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where K is the Wigner or time- inversion operator. In Refs. [|5] and || it is proved that 
this operator changes z into z* so really a* can be considered as a function of z* rather 
than z. 

Now if we make a Bogolyubov transformation among operators labelled by z, the 
cannonical conmutation relations (3.6) must be kept invariant. Therefore the usual 
relations among Bogolyubov coefficients, such that = aiSij, = /3i5ij, namely 



in the complex case it reads 



M 2 -|A/ = 1, (B.4) 



a z a* z ,-(3 z (3:, = l. (B.5) 



Now let us study the functions A 2 = (3 z /a z . The poles of these functions originate 
the poles of the S-matrix. Let us compute: 

\mA ici 



\a z aU \l+p g /3*.\ |1 + C| 



— < 1, (B.6) 

where P z f3*, = ( = £ + irj. Then function f3 z f3** / 'a z a*, has a bounded modulus and 
therefore it cannot be analytic in all the complex plane if it is not a constant. But a 
constant is not a generic function. Therefore function A 2 A 2 , must have some singularities 
some poles and/or some cuts. If it is a cut, either A z or A* has a cut and therefore also 
the S-matrix, but a cut can be considered (in all the demonstrations of this paper) like 
an infinite set of poles and the process is finished and the theorem is demostrated. On 
the other hand let us suppose that the function A 2 A 2 , has just a pole at z$. Therefore 
either A 2 has a pole at z or A*, has a pole at this point. Thus either S has a pole at Zq 
or at Zq, thus in any case S has at least one pole. 

Now as function A z has a pole let us say at z , we can define a function 

A« = A Z - \ ResA ^ . (B.7) 

Z- Zq 

Let us suppose that this function has no poles. We can compute the modulus 
|A^A^*| which obviously is bounded since 

|A«A^| < 1, (B.8) 

if z and z* are far enough of z , and also near to z since AW has no pole there, because 
the only pole of A 2 has at zq and it has been eliminated via the second term of the r.h.s. 
of (B.7). Then function A^A^t is analytic and bounded in all the complex plane and 
therefore it is a constant. But again a constant is not a generic function so A^A^* 
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have some singularities, either a cut, in which case the proof is finished, or it has a pole 



at Zq~\ Now we can define a new function 



A« = A?> - (B.9) 



and so forth and then prove that A 2 , and therefore S, have both an infinite set of 
poles, or a cut, which can also be considered as an infinite continuous set of point-like 
singularities, for all the issues considered in this paper. We can also introduce another 
way to forsee the presence of this infinite set of poles. 

According to what we have just said the equation a ik = of (3.4) has at least a root 
(since we cannot say that (3 ik would have any pole in a generic case). But this equation 
has a "periodical" nature (in the complex plane) as we can see if we try to solve a similar 
equation like 

= £• (B.10) 



r(i-s) 

Using (1.2.2) of Ref. [16] this equation can be written as 

-T(z)sm7rz = f. (B.ll) 

7T 

Let us only consider the case £ = yE 1Z, z = x<E 1Z, x > then 

-r(x)simrx = y. (B.12) 

7T 

The function T(x) is always > if x > and it is a growing function. Then the 
l.h.s. of (B.12) is a periodical function modulated by a growing one, then we can see the 
periodical nature of the problem. Then the function (B.12) intersects at an infinite set 
of points the straight line y = const. Therefore (B.12) has an infinite number of roots. 
As oiik = is a generalization of (B.12) and it has at least a roots, we can conclude that 
it has an infinite number of root owing its periodical nature. 
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